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ABSTRACT 

Very high-energy (VHE; E>100 GeV) y-rays have been detected from a wide range of astronomical objects, such as pulsar wind 
nebulae (PWNe), supernova remnants (SNRs), giant molecular clouds, y-ray binaries, the Galactic Center, active galactic nuclei 
(AGN), radio galaxies, starburst galaxies, and possibly star-forming regions as well. At lower energies, observations using the Large 
Area Telescope (LAT) onboard Fermi provide a rich set of data which can be used to study the behavior of cosmic accelerators in the 
MeV to TeV energy bands. In particular, the improved angular resolution of current telescopes in both bands compared to previous 
instruments significantly reduces source confusion and facilitates the identification of associated counterparts at lower energies. In 
this paper, a comprehensive search for VHE y-ray sources which are spatially coincident with Galactic FermifLKT bright sources 
is performed, and the available MeV to TeV spectra of coincident sources are compared. It is found that bright LAT GeV sources 
are correlated with TeV sources, in contrast to previous studies using EGRET data. Moreover, a single spectral component seems 
unable to describe the MeV to TeV spectra of many coincident GeV/TeV sources. It has been suggested that y-ray pulsars may be 
accompanied by VHE y-ray emitting nebulae, a hypothesis that can be tested with VHE observations of these pulsars. 

Key words. Gamma rays: observations - Galaxy: general - (Stars) pulsars: general - ISM: supernova remnants - X-rays: binaries 



1. Introduction 

Our understanding of the very high-energy (VHE; E >100 GeV) sky 
has greatly improved during the last few years, thanks to the high sen- 
sitivity of current imaging atmospheric Cherenkov telescopes (IACTs), 
e.g. H.E.S.S., MAGIC, and VERITAS. They typically cover the energy 
range of ~100 GeV up to several tens of TeV, and provide an angular 
resolution of ~ 6'. This allows spectral and morphological studies of 
the various types of VHE sources: pulsar wind nebulae (PWNe), su- 
pernova remnants (SNRs), giant molecular clouds, y-ray binaries, the 
Galactic Center, active galactic nuclei (AGN), radio g alaxies, starburst 
galaxie s, and possibly star-forming regions as well. See Aha ronian et alJ 
( 2008eJ) for a revi ew of the field i n 2008 , with more recent updat e s give n 
by the H.E.S.S. dChaves et all l2009tJ) MAGIC dZanin et all [2009). 
and VERITAS collaborations dOng et alll2009tlWeinstein et aU r2009). 
However, many of the sources are not yet identified at other wave- 
lengths, e.g., nearly a third of the Galactic H.E.S.S. sources have no firm 
identification; in many cases, there are multiple plausible counterparts 
while, in other cases, no viable counterparts have yet been identified. 

Gamma-ray observations of Galactic sources can help us to solve a 
number of important astrophysical questions, including (1) the physics 
of pulsars, PWN and SNR; and (2) the origin of cosmic rays. Our 
Galaxy contains a large number of cosmic accelerators, where particles 
are accelerated to highly-relativistic energies (up to at least 10 14 eV). 
The origin of cosmic rays is still not well known, largely because of 
the lack of directional information of these particles. These very en- 
ergetic particles can be traced within our Galaxy by a combination of 
non-thermal X-ray emission and y-ray emission via leptonic (such as in- 
verse Compton scattering of electrons, Bremsstrahlung and synchrotron 
radiation) or hadronic (via the decay of charged and neutral pions, due 
to interactions of energetic hadrons) processes. Therefore, observations 



of y-rays at energies ^100 MeV can probe the sources of particle accel- 
eration. 

The Large Area Telescope (LAT), onboard the Fermi Gamma-ray 
Space Telescope, provides the best information of the non-thermal sky 
in the energy range from 20 MeV to 300 GeV. The point-source sen- 
sitivity of LAT is ~10~ 8 p hcrrr 2 s -1 above 100 MeV in one year of 
survey-mode observations (Atw ood et all 12009). which is an order of 
magnitude better than that of its predecessor, the Energetic Gamma 
Ray Experiment Telescope (EGRET). Its angular resolution is 5j 0:6 
above 1 GeV, which is particularly important for identifying y-ray 
sourc es with multi-wave length counterparts and revealing their na- 
ture ( I Atwood eta l., 2009). As an important step towards the first source 
catalog, the LAT collaboration has published a bright source list (BSL) 
that includes 205 sources, designated with the pre fix OFGL, using dat a 
taken during the first three months of observations dAbdo et al., 2009a). 
Among them, 121 sources are identified with AGN and one with the 
Large Magellanic Cloud. Most of the remaining 83 sources are believed 
to originate in our own Galaxy. It is natural to investigate which of them 
also have been detected at energies ^100 GeV. 

The search for VHE counterparts of LAT sources is important for 
the following reasons: 



it aids the identification of the true nature of the LAT sources 
through their VHE counterparts; 

for pulsars, it helps us to identify their VHE-emitting nebulae; 
it may provide us with broad-band y-ray spectra, thereby better 
constraining the emission mechanisms (e.g. distinguish between 
hadronic and leptonic scenarios). 



* now at Institute of Astronomy and Department of Physics, National 
Tsing Hua University, Hsinchu 30013, Taiwan 



iFunk et al.1 (2008) compar e y-ray sources in the third EGRET 
(3EG) catalog faartman etafl. 1 1 9990 and the 22 H.E.S.S. sources 
kn own at the time within th e region of / = -30° to 30°, b = -3° to 
3° dAharonian et all 1200611) . They do not find spatial correlation be- 
tween the two populations. Though some coincidence cases are found, 
the authors conclude that these few cases can be explained by chance 
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coincidence. However, due to the capabilities of EGRET, this study suf- 
fers from the following limitations: (1) The sensitivity of EGRET is 
lower than that of LAT. The lack of photon statistics leads to poorly 
constrained spectral indices and the spectra terminate S10 GeV at the 
upper end for a typical source; (2) EGRET sources are only localized 
at degree scales, which is much larger than the angular resolution of 
IACTs. The second point is the instrumental reason which explains th e 
weak correlation of EGRET and H.E.S.S. sources dFunket aUl2008l) . 
These shortcomings are now largely overcome due to the enhanced per- 
formance of LAT over EGRET. In addition to the above caveats, they 
do not consider the extension of the VHE y-ray sources in their anal- 
ysis. As such, the full potential of this search is not realized for very 
extended sources lik e the SNR RX J 1 7 13. 7-3946, as pointed out by 
Tiboll a et all d2009a ). After the launch of LAT, one largely benefits 
from the increased LAT an gular resolution over previous studies. As 
noted in Atwood et al. (2009), EGRET could not distinguish the GeV 
emission of RX J1713.7-3946 from 3EG J1714-3857, while the ca- 
pabilities of LAT allow one to study individual sourc es in this region, 
which contains three VHE y-ray sources (See Fig. 1 in Aharonian et al., 
I2008bh . 

The water Cherenkov detector MILAGRO covers the energies 
above ~1 TeV and its angular resolution can reach < 1°. Using 
MILAGRO, a search for y-rays f rom the Galactic LAT BSL has been 
performed bv lAbdo et"aU d2009g|) . They found that 14 sources (of the 
selected 34) show evidence for multi-TeV y-ray emission at a signifi- 
cance of > 3tr, although most of the source can didates cannot be estab- 
lished as firm detection on an individual basis dAbdo et all 120093) . 

In this paper, a search for VHE co unte rparts of all the pr esumed 
Galactic sources in Abd oltal] d2009ah and|Abdo et al. C 2010h is per- 
formed, with spatial coincidence as the primary criterium for associa- 
tion. The extensions of the VHE y-ray sources are taken into account 
and the search is not limited to the H.E.S.S. Galactic plane survey re- 
gion. The broad-band MeV to TeV spectra of coincident sources are 
then presented. 

2. Search of spatial coincidence 

2.1. The Fermi and VHE catalogs 

Abdo et al. (2009a) present 205 point-like sources which were detected 
at or above the 10-cr level in the 0.2-100 GeV band, based on three 
months of observations (August 4, 2008 - October 30, 2008). The list 
is not flux-limited and thus is not uniform. The following information 
is given for each source: its position, positional uncertainty (95% confi- 
dence level, C.L.), significance, flux in two energy bands (100 MeV- 
1 GeV and 1-100 GeV), and whether or not there is e vidence for 
variabil ity over the above-mentioned period. In addition, Abdo et al. 
(2009a) assign for each source its source class, as well as y-ray and 
lower energy association (if any). Those sources that are classified as 
extragalactic (all AGN and the Large Magellanic Cloud) are not con- 
sidered in this work. 

The remaining source list contains 83 sources, comprising 15 
radio/X-ray pulsars, 15 pulsars newly discovered by the LAT, two high- 
mass X-ray binaries (HMXBs), one globular cluster (47 Tucanae), 13 
SNR/PWN candidate^ and 37 sources without ob vious counterparts a t 
lower energies (a mong them the Galactic Center: lAbdo et al.. 2009a). 
Abdo et alj OOlOh presents the first LAT pulsar catalog. Those 16 pul- 
sars that are not present in the above BSL are also included in this 
study. Therefore, most of the LAT bright sources considered in this 
work should be Galactic in origin. 

T he number of VHE y-ray sources is larger than 50 as of Fal l 
2009 dAharonian et al.Ll2008d:IChaves et alll2009al : IZanin et alll2009l : 
lOng et aU?l2^0^nWem^teTn^t^lTl2"009h . Galactic sources summarized 
in the above publications are used in the coincidence search in this 
work. Therefore, our comparison is based on published sources only. 
The VHE y-ray source positions and source extensions in this work 
are taken from the corresponding publications shown in Tables [T] [2] 

1 possibly associated with SNR or PWN, but the emission may come 
from unknown pulsars 



and [3] At higher energies, the MILAGRO collaboration reported ev- 
idenc e for multi-TeV y-r ay emission from several LAT source posi- 
tions (Abdo et al., 2009a). Only those source candidates with a signifi- 
cance larger than 5tr are regarded as TeV sources here and are included 
in this studjQ. With several tens of known sources in both the GeV and 
TeV bands, a systematic cross-correlation study is conducted. 

2.2. Level of spatial coincidence 

To quantify the level of spatial coincidence, the following scheme is 
employed. Let d be the distance between a centroid position best fit by 
LAT and the best-fit centroid of a nearby VHE y-ray source. The radius 
of the 95% confidence region for the LAT source is the uncertainty on 
the centroid position of the given LAT source, which is typically ~ (FA. 
Most VHE y-ray sources are extended, with a typical size of 0:1 - 0°5. 
Let e be the sum of the radius of the 95% confidence region and the size 
of the VHE y-ray source. 

For each LAT source, if a VHE y-ray source was found where 
d - e < 0, the source pair is classified as a spatially coincident case (i.e. 
category Y - Yes). Given the possible extended nature of many LAT 
bright sources, a category P (for Possible) is defined for pairs where 
< d - e < 0.3°, so that the sources seen by LAT and the VHE in- 
struments may actually overlap each other; they are possibly coincident 
cases. If no reported VHE y-ray source was found with d - e < 0.3°, 
the LAT source falls into the coincidence level N (for No), i.e., no co- 
incidence with any VHE y-ray source. If there are several nearby VHE 
y-ray sources, only the closest VHE y-ray source would be considered. 

2.3. Spatial coincidence GeV/TeV pairs 

In the search, 24 coincident cases (Y, among them two are coincident 
with MILAGRO source only) and 7 possibly-coincident cases (P) are 
found. The results are presented in Tables[T][2] and[3] No reported VHE 
y-ray sources are found in the remaining 68 sources. 

According to the nature of the LAT sources, the results are summa- 
rized as follows: 

1. Eight LAT pulsars are spatially coincident with a source detected 
using IACTs, which may be the VHE-emitting PWN. There are 
two additional pulsars which are possibly coincident with an IACT 
source. Two others have a MILAGRO counterpart, but have not yet 
been detected by IACTs. 

2. Among the 13 SNR/PWN candidates in the Fermi BSL, more than 
a half (7) are spatially coincident with a VHE y-ray source, and 
another one is a possibly coincident case. The seemingly high frac- 
tion of coincidence is partly due to a better coverage of the inner 
Galaxy region, where most SNR/PWN candidates are located. This 
results in a generally better sensitivity for this class of sources than 
for other classes. 

3. The two HMXBs listed in the BSL 
(0FGL 10240.3+61 13/LS I +61 303 and 
0FGL J1826.3-1451/LS 5039) are both found to be spatially 
coincident with a VHE y-ray source, identified with the same 
corresponding binary. 

4. Five of the 37 unidentified 0FGL sources are spatially coincident 
with a VHE y-ray source. The number increases to nine if possibly 
coincident cases are included. 

In addition, we note that a new VHE source near PSR Jl 1 19-6127 
has been announced in an oral presentatior0, but it has not been pub- 
lished with any written reference to our knowledge. Given that the best- 
fit centroid and extension were not given, we do not treat it as an entry 
in our samplfl 

2 For example, HESS 11833-105 dDiannati-Atai et all 120 08b). hav- 
ing a significance of 4.0cr only but included in dAbdo et alll2oToh as a 
counterpart of LAT pulsar PSR J1833-1034, is not considered here. 

3 See http://cxc.harvard.edu/cdo/snr09/pres/DjannatiAtai_Arache_v2.pdf 

4 In the first LAT catalog, that can be found at 
http://fermi.gsfc.nasa.gov/ssc/data/access/lat/lyr_catalog/, the authors 
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With such a large number of coincident cases, the relationship be- 
tween the GeV and TeV sources is explored. In the next section, the 
y-ray spectral energy distributions (SEDs) are constructed for those co- 
incident and possibly coincident GeV/TeV source pairs with published 
VHE spectrum. 



3. The y-ray spectral energy distributions 

3.1. Construction of PL spectrum in the LAT energy range 

lAbdo et ail d2009al) provide photon fluxes and respective errors in 
two energy bands: low energy (100 MeV-1 GeV) and high energy 
(1 - 100 GeV). Since photon spectra are not given in the BSL, we at- 
tempt to estimate the spectra of individual sources. 

Assuming that a single pure power law (PL) represents the spectrum 
in the LAT energy range, the photon flux in the low (10 2 -10 3 MeV) 
and high energy (10 3 -10 5 MeV) bands, respectively, are given by 

F23 = k j Q1 E~ T dE and F35 = k E~ T dE, where E is measured 
in GeV, T is the photon index, and k is the normalization at 1 GeV. 
From these two expressions, k and T can be calculated. Using the avail- 
able flux errors (AF23 and AF 35 ), uncertainties in k and T (Ak and Ar) 
are obtained by error propagation. The spectra are then constructed in 
the form of "bowties". For those where F23 is given as a 2-cr upper 
limit, the calculated T can be treated as an upper limit and the recon- 
structed spectra can be seen as the "softest possible" PL spectra. The 
PL spectra are plotted from 100 MeV up to a certain maximum en- 
ergy, F max (< 100 GeV), which is defined by requiring that the photon 
spectrum above F max contains 10 photons over the three months of ob- 
servationfl This results in a range of values for F max from ~3 GeV to 
100 GeV. The single PL assumption does not hold in general. Given the 
limited information available in the BSL, such an assumption should 
be seen as a very rough estimation of the source spectra and it is 
used in this work for a visual GeV/TeV spectral comparison. A cut- 
off between the GeV and TeV bands has been measured for several 
sources including pulsars. Therefore, we also plot the best-fit spec- 
tra when a detailed LAT spectrum is available in the literature (Vela, 
Crab, Geminga, PSR J1706-44, and LS I +61 303). For the cases of 
0FGL J0617.4+2234 and 0FGL J1746.0-2900, the double PL spec- 
tra derived for 3EG J0 617+2238 and 3EG J1746-2851, respectively, 
bv lBertsch e t al. (2000) are also shown for comparison. 

3.2. The MeV-TeV SEDs 

The sources considered here are those 0FGL/VHE pairs with spatial 
coincidence levels Y and P and with VHE spectral information available 
in the literature. In the case of HESS J1923+141 where only a VHE flux 
is given, a typical spectral index is assumed. In addition, there are two 
pulsars for which a MILAGRO candidate counterpart is reported but 
there is no VHE y-ray detection using IACTs (see Table[2). 

SEDs of the 28 cases in the energy range from 100 MeV to 
>1 TeV are depicted in Figs. [T] to [9] Systematic errors in spectral in- 
dices and normalization are not shown, which for TeV spectra are 
~20% for most sources and for GeV spectra ar e ~ 20-30% (the lat - 
ter is inferred from flux estimation systematics in lAbdo et all 12009a). 
An overall inspection of the SEDs immediately shows that single PLs 
from 200 MeV to ~10 TeV cannot describe most GeV-TeV y-ray 
spectra. This is not surprising given the large range in photon en- 
ergy (i.e. five orders of magnitude), as no photon spectrum from any 
emission mechanism is expected to be unbroken for such a large en- 
ergy span. The only example for which a pure PL may still work is 
0FGL J1836.1-0727/HESS J1837-069, which is a possibly coincident 
pair (P). The most common board-band behaviors are a cut-off at ener- 
gies below ~100 GeV (dominating in the pulsar class), and a spectral 



claim that there is an association of the LAT source 1FGL J 1 1 19.4- 
06127 with a VHE source, designated HESS Jl 119-614, which may 
be the same VHE source. 

5 using the LAT on-axis effective ar ea above 1 GeV of ~8000 cm 2 
and a mean on-axis exposure of ~1 Ms dAbdo et a l.. 2009a) 



break between the LAT and the VHE bands (dominating in the uniden- 
tified LAT sources). 

The SEDs of the LAT source classes including pulsars, SNR/PWN 
candidates, and unidentified y-ray sources are presented in this section. 
LS I +61 303 and the Galactic Center region are discussed in Sect. 14.51 
and !4.6l respectively. 

3.2.1. Pulsars 

Figure [T] shows the four y-ray pulsars known in the EGRET era, Fig. [2] 
shows the four radio pulsars first detected in y-rays by LAT, and Fig. [3] 
shows the three new pulsars after a blind search for pulsations in the 
LAT data (Abdo et al., 2009b). Besides the Crab, no off-pulse emission 
is found in the LAT data of the other 10 pulsars, suggesting that most of 
the emission from pulsars seen with LAT is pulsed and originates from 
the pulsars themselves. On the other hand, extended regions are seen at 
energies above 100 GeV in these 10 cases (except for the Crab, which 
appears point-like to all IACTs). Their VHE emission (>100 GeV) is 
unpulsed and for many of them (e.g. Vela X) this emission have been 
attributed to PWNe, although in some cases there exists other possibil- 
ities to explain the VHE y-ray source (e.g. a spatially coincident SNR). 

The SEDs of the pulsars essentially depict the pulsed component in 
the LAT energy band and the unpulsed component in the VHE band. 
It is suggested that the emission below and above ~100 GeV comes 
from two different emission regions, e.g. pulsed emission from the pul- 
sar magnetosphere and unpulsed emission from the PWN. It can be seen 
that (1) a cut-off must exist between the LAT "bowties" and the corre- 
sponding VHE spectra (with the notable exception of the C rab) - this is 
demo nstrated with a detailed spectral study of pulsars (e.g., Abd o et all 
2010); (2) the energy output at GeV energies is at least an order of 
magnitude higher than that in the VHE band. This indicates that for 
the pulsar population presented in this section, the PWNe radiate less 
energy than the y-ray pulse emitting regions. 

However, the PL derived LAT spectra are not always good repre- 
sentations of the reported spectra for individual sources. This is demon- 
st rated in Fig.[pwher e both the "bowtie" spectra and the derived spectra 
in I Abdo et al] d2010h are shown. In all the other cases, only the spectra 
as presented in lAbdo et al] (2010) are depicted. 

3.2.2. SNR/PWN candidates 

The SEDs of those 0FGL sources classified as SNR/PWN candidates 
are shown in Figs. [4] and [5] The GeV-TeV spectral connection varies 
among the sources in this class. The TeV spectra are not simply the 
PL tails of the GeV spectra. There are cases where the extrapolation 
of the LAT "bowtie" to TeV energies is at least an order of mag- 
nitude higher than the measured VHE flux level (e.g. the spatially 
coincident case 0FGL J1801.6-2327/HESS J1801-233, a cut-off oc- 
curs between the two energy bands), while for another coincident case 
(0FGL J1834.4-0841/HESS J1834-087) the PL extrapolation to the 
VHE band is below the measured VHE level and a second spectral com- 
ponent above ~200 GeV is needed to explain the TeV excess. 

There is only one case (0FGL J0617.4+2234) where a bro- 
ken PL describes the LAT spectrum better than a single power 
law. The "bowties", which are derived a priori from PLs, may 
be closer to the real spectra compared to the case of pulsars. 
If that is the case for 0FGL J1801.6-2327/HESS J1801-233 and 
0FGL J1923.0+1411/HESS J1923+141, a spectral break may occur at 
energies in the largely unexplored energy range of 10-100 GeV for 
these two sourcefl 

3.2.3. Unidentified LAT sources 

The SEDs of those 0FGL sources without obvious counterparts are 
shown in Figs. [6] and [7] There is so far no published spectra of this LAT 



6 The LAT spectrum for 0FGL J 180 1.6-2327 is the softest possible 
PL, while the HESS J1923+141 spectrum is derived assuming a PL 
index T = 2.8. 
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source class. For the case of OFGL J1839.0-0549/HESS J1841-055, 
the spectrum may span from ~100 MeV to ~80 TeV, with a pos- 
sible break within or close to the "energy gap" at ~60-500 GeV If 
the GeV and TeV sources are indeed associated, they might repre- 
sent a group of "dark accelerators" which have a broad y-ray spec- 
trum. All SEDs are consistent with the assumption that a spectral 
break exists between the two energy bands, except for the case of 
OFGL J1805.3-2138/HESS J1804-216, a spatially coincident case (Y). 

3.3. Comparison of the flux and photon indices in the GeV 
and TeV energy bands 

A comparison of the flux levels in the GeV and TeV energy bands for 
coincident GeV/TeV sources (category Y) is depicted. Figure [Tolshows 
the photon flux in the 100 MeV - 1 GeV band plotted against that in 
the 1-10 TeV band (derived according to Sect. ). For most sources, the 
photon flux in the 1-10 TeV band, F[_ 10 tcV> is about 10~ 5 to 10 -6 that in 
the 0.1-1 GeV band. Figure[TT]shows the photon flux in the 1-100 GeV 
band plotted against that in the 1-10 TeV band. For most sources, pho- 
ton flux in the 1-10 TeV band, Fi_iotcV, is about 10~ 4 to 10~ 5 that in 
the 1-100 GeV band. 

Figure [T2l depicts the photon indices in the 0.1-100 GeV band de- 
rived according to Sect. l3.3l against photon index in the 1-10 TeV band. 
It can be seen that the TeV spectra are similar to or harder than the GeV 
spectra for most sources, i.e. < (Fi_ 10 Tev _ r .i_iooGev) L 



4. Notes on selected sources 

Although detailed analysis of the LAT data for each individual source 
is beyond the scope of this paper, some comments on the following 
sources are given. 

4. 1 . Crab pulsar and nebula 

The Crab pulsar and nebula are among the best-studied non-thermal 
objects in the sky. The pulsed emission above 100 MeV and up to 
~10 GeV is clearly detected with LAT. Two strong peaks are seen in the 
phase histogram. A spectral fit of the pulsed emission usi ng a PL with an 
expon ential cut-off gives a cut-off energy of ~8.8 GeV iGrondi n et all 
2009). There is evidence of pulsed emission u p to ~25 GeV, as mea- 
sured using the MAGIC telescope (see Fig. [Tl lAliu et all l2008h . The 
flux reported by MAGIC is consistent with the exponential cut-off in 
the spectrum measured by LAT. 

Evidence for unpulsed emission was already present in the EGRET 
data dde Jager et all [1996). The LAT measurement of this component 
can be well fit by a single PL with T ~ 1.9 up to ~300 GeV. This 
unpulsed spectrum agrees well with the VHE spectra measure d by the 
IACTs MAGIC, H.E.S.S. and VERITAS dGrondin et aill2009h . In par- 
ticular, there appears to be a deviation from a pure PL in the MAGIC 
spectrum below -100 GeV dAlbert et alll2008l) . 

4.2. Vela pulsar and Vela X 

The Vela pulsar is the strongest persistent GeV source and was the first 
target of LAT observations. The complex pulse profile is dominated 
by two peaks with a pronounced "bridge" between them. The phase- 
averaged spectrum, which is essentially the pulsed emission, can be 
well described by a PL with an exponential cut-off at ~2.9 GeV. The 
off-pulse emission is much weaker, and a 95% C.L. upper limit of the 
photon flux of 1.8 x 10~ 7 cirT 2 s is derived at the puls ar position in the 
0.1-10 GeV band (shown in Fig.fTl lAbdo etaill2009c[) . 

To the south of the pulsar, an extended VHE y-ray source spatially 
coincident with the Vela X region, HESS J0835-455, has been detected. 
The observations r epresent the first measur ement of a SED peak in a 
VHE y-ray source ( Aharo nian et all 12006b). The PL with exponential 
cut-off fit of this PWN is reproduced in Fig.[TJ An analysis of the Vela X 
region does not establish a n ebula component based on the first three 
months of LAT observations ( Lemoine -Goumard et alll2009h . 



4.3. Geminga 

The Geminga pulsar is t he first known radio-quiet y-ray pulsar in the 
sky dBertsch et allll992l) . See Fig. [T] for its SED. While EGRET data 
are well fit by a single PL up to 2 GeV (but show s evidence of a cut- 
off above 2 GeV; Maver-Hasselwande r et all 1 19941) . the cut-off energy 
is de termined to be ~2 .6 GeV using the first seven months of LAT 
data dCeliket all 120091) . There appears to be an excess at ~20 GeV 
when compared to the fit with a PL with exponential cut-off. The reason 
may be due to low statistics or effects of the fitting method, but it might 
also indicate a seperate and harder spectral component (Celiketal., 
2009). There is as yet no evidence for unpulsed emission. 

Evidence for multi-TeV emission around the pulsar was rep orted in 
the MILAGRO survey of the Galactic plane dAbdo et al l 120071) and in 
the sea rch for MILAGRO counterparts of Fermi sources dAbdo et all 
2009g), using a point source analysis, at ~3<x (post-trial) significance 
levels. Assuming that the emission is extended, the significance in- 
creases to 6.3<t at the position of the pulsar. If the detection is real, the 
size of the MILAGRO emission is ~2°6. At a distance of only ~ 250 pc, 
this extent is similar to more distant PWN dAbdo et all l2009g|) . On the 
other hand, VERITAS observations resulted in no detection but rather a 
99% C.L. flux upper limit (above 300 GeV) of 2 x 10' 12 cm' 2 s" 1 , as - 
suming point source emission from the pulsar (Finnegan et al., 2009). 
Although IACTs suffer from reduced sensitivity when observing very 
extended source (which scales as 6T 1 with 8 being the extension), obser- 
vations of Geminga with IACTs are crucial for verifying the MILAGRO 
claim and helping us to understand the y-ray emission from Geminga. 

4.4. PSRB1 706-44 

Gamma-ray pulsations from PSR B 1706-44 were discovered by 
EGRET; the observations revealed a triple-peaked pulse profile but 
no evidence for unpulsed emission 1 Tho mpson et all 1996). More re- 
cently, the pulsar was also detected by FermifLKI as the bright source 
OFGL J1709. 7-4428. The phase-averaged spectra measured by EGRET 
and LAT are both well-described by a broken PL (up to 30 GeV, 
in the case of the LAT spectrum). The break energy measured by 
LAT is 3 GeV, while in deriving the EGRET sp ectrum, it is fixed at 
1 GeV dThompson et allll996l ; rBertsch et alll2000h . The LAT PL index 
steepens from a value of ~ 1.9 (below 3 GeV) to ~3.3 (above 3 GeV), as 
shown in Fig. [TJ This spectrum and the PL spectrum derived using the 
method described in Sec. l3.3l are both consis t ent wit h the photon flux in 
the 1-100 GeV band reported in lAbdo et all d2009al) . 

The discovery of an extended source of VHE emission in the vicin- 
ity of PSR B 1706-44 was recently reported by H.E.S.S. dHo ppe et all 
2009). The TeV source is quite hard (r ~ 2.0), more so than the high- 
energy part of the pulsar spectrum. The VHE y-ray emission is sug- 
gested to be relate d either to a relic PW N of PSR B 1706-44 and/or the 
SNR G343.1-2.3 dHoppe et alll2009l) . 

4.5. LS I +61 303 

LS I +61 303 is the first X-ray b inary where periodi c y-ray emission 
has been detected at both GeV (|Abdo et all l2009d) and TeV ener- 
gies dAcciari et ail 12009 j ; lAlbert et all 120091) . Its SED is shown in 
Fig. [8] The "bowtie" looks nicely connected to the measured VHE spec - 
trum, but a cut-off energy at ~6 GeV is reported dAbdo et a l.. 2009d). 
Furthermore, the timing measurements in both the GeV and TeV bands 
show that the maximum emission occurs at different orbital phases, 
namely, close to periastron for <100 GeV emission and close to apas- 
tron for VHE emission. Thi s suggests different emission mechanisms in 
the two bands, as noted in dAbdo e t al., 2009d). 

4.6. Galactic Center Region 

The Galactic Center is among the richest and most complex regions in 
the Galaxy, due to the large number of possible sources and the dif- 
ficulty of correctly modeling the diffuse emission due to cosmic-ray 
interaction with the local molecular clouds. This problem is extremely 
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relevant at GeV energies, as demonstrated by EGRET measurements. 
The discovery of new VHE y-ray sources close to the Galactic Center 
is important for studying the role of diffuse Gal actic emission vers us the 
emission from resolved sources in this region dTibolla et alil2009l) . 

One GeV source, OFGL J 1746.0-2900, is detected with a sig- 
nificance of 36cr in the neighborhood of the Galactic Center. The 
best-fit position for OFGL J1746.0-2900 is R.A.= 17 h 46 m l'.'4, Decl.= 
-29°0 , 18" (J2000) with a 95% C.L. error radius of 4' dAbdo et all 
2009a). The H.E.S.S. Collaboration also reported a detection of a 
source towards the Galactic Center, localized at R.A.= 17 h 45 m 39'.'6 ± 
q'4(stat) ± g'4(sy s), Decl.= -29°0'22" ± 6" (stat) ± 6" (sys) (J2000; 
lAcero et"alT I2010T) . Based on the procedure described in Sect. 12.21 the 
OFGL J1746.0-2900/HESS J1745-290 pair falls into the category of 
possibly coincident cases. With better photon statistics, one of the fun- 
damental questions that the LAT can hopefully address is whether or 
not the GeV and TeV sources are indeed spatially conincident. 

The spectra of OFGL J1746.0-2900 and HESS J1745-290 are 
shown in Fig. [9] The spectra in the two bands do not appear to be well 
described by a single PL, and there seems to be an order-of-magnitude 
drop-off in the energy range ~10-100GeV. Although detailed analysis 
of the LAT data is out of the scope of this paper, this simple inspec- 
tion does not indicate that they are from the same emission component 
(although large uncertainties due to systematics in this region do not 
permit stronger conclus ions at this time). Fo r reference, the broken PL 
fit of 3EG J1746-2851 dBertsch et alll200Cl) is also shown in Fig.[9] 

5. Discussion 

In this work, the first comparison of the GeV and VHE y-ray sources 
after the launch of LAT is presented, which takes the advantage of the 
significantly improved LAT angular resolution and sensitivity compared 
to EGRET. 

Below are a list of preliminary results drawn from this work: 

1. With the better localization and morphological information of 
VHE y-ray sources compared with OFGL sources, the nature 
of the LAT sources may be better revealed through their VHE 
counterparts. Table [T] lists the potential counterparts of some 
VHE y-ray sources which are coincident with OFGL sources. 
For example, HESS J1804-216 may be related to W 30, which 
may help in understanding the nature of the unidentified source 
OFGL J1805. 3-2138. 

2. Results of several LAT-detected pulsars sh ow cut-offs at ene rgies 
~ 1—10 GeV, similar to the assessment of iFunk et al.l d2008l) for 
EGRET-detected pulsar systems. Therefore, a VHE counterpart 
(~0.1-10 TeV) of a LAT pulsar most likely represents the asso- 
ciated PWN, with a shell-type supernova being a viable alterna- 
tive. This is particularly important for those new pulsars discovered 
by LAT. The VHE counterparts coincident with the six LAT pul- 
sars may be the associated PWN, although other explanations (e.g. 
shell-type SNR) are also possible. The question of whether or not 
typical y-ray pulsars are accompanied by VHE-emitting nebulae 
can be tested by observing them in the VHE domain. 

3. Through broad-band y-ray spectra of SNRs, one may in princi- 
ple distinguish between hadronic and leptonic scenarios. A study 
of RX J1713. 7-3946 usin g five years of simulated LAT observa- 
tions dAtwood et all [2009) shows that the energy flux level for the 
hadronic scenario differs by around a factor of two from that for 
the leptonic scenario and that a spectral break may be more promi- 
nent for the latter. The SNR sample shown in Figs. [4] and [5] do not 
seem to support either scenario, although it is too early to draw any 
conclusion based on the three-month OFGL dataset. If a hadronic 
scenario is found to be more viable, this would support the current 
hypothesis that shell-type SNRs are cosmic-ray sources. 

4. Previous studies did not reveal a strong co rrelation between 
the GeV/TeV populations. iReimer et al.1 (2008) list 16 H.E.S.S. 
sources without counterparts from the 3EG catalog. Among 
them, new associations are found in the present study and 
are presented in Table [4] tha nks mostly t o the discovery of 
new GeV sources with LAT. IReimer et al.l d2008h also present 



11 sources in the 3EG catalog without H.E.S.S. counterpart. 
Among them, OF GL J1709.7-4428 ( the OFGL counterpart of 
3EG J1710-4439; lAbdo et all l2009al) is now found to be as- 
s ociated with HESS J1708-443, a source discovery reported 
in lHoppe et alJ J200l^ . 

5. All spatially coincident GeV and TeV pairs during the EGRET era 
are essentia lly consistent wit h one single spectral component (See 
Figs. 4-6 in Funk et al., 2008). With the significantly enhanced sen- 
sitivity of LAT, new relations between the GeV and TeV spectra 
are apparent in the SEDs. The SNR candidate OFGL J1834.4-0841 
and the unidentified OFGL J1805.3-2138 (and their likely VHE 
counterparts) represent the first examples for which the GeV/TeV 
spectrum cannot be treated as a single emis s ion com ponent. A sim- 
ilar conclusion is reached by lAbdo et aU d2009dl) for a HMXB 
(LS 1 +61 303), based on the light curves and spectral incompat- 
ibility of th i s source in the two bands. 

6. lAbdo et ail d2009d) consider a probability that many unidentified 
LAT sources are extragalactic, so as to explain the low rate of find- 
ing coincident MILAGRO emission among the unidentified LAT 
sources. This idea might also explain the non-detection of VHE 
counterparts of a majority of the unidentified LAT sources. On the 
other hand, the extended nature of all the five spatially coincident 
cases (HESS J1023-575, HESS J1804-216, HESS J1841-055, 
HESS J1843-033, HESS J1848-018; if proved to be real associ- 
ation) would exclude an extragalactic origin of the corresponding 
LAT sources. 

7. Although VHE observations only cover a small part of the whole 
sky, they do cover the majority of the inner Galaxy, e.g., the 
H.E.S.S. telesco pes have surveyed th e region of I = -85° to 60°, 
b = -3° to 3° (Chav eTet all 120 09a). In this region, there are 41 
Fermi bright sources. Among them, 16 are found coincident with 
a VHE counterpart. This fraction (~2/5) is higher than that for 
EGRET where about 1/4 of the EGRET sources (in a smaller re- 
gion) are found to have a coincident VHE counterpart ( Fun k et all 
2008). Moreover, the number rises to 21 (out of 41) if possibly co- 
incident cases are included and the fraction becomes 50%. The LAT 
radii of the 95% confidence region are in general much smaller than 
the EGRET error boxes, which further strengthens the case of a 
higher fraction for LAT. Even though the VHE extensio n is taken 
into account in this study (but not in Funk et al., 2008), a typical 
extension is of the same order as a LAT positional uncertainty. A 
breakdown of the number of coincidence cases for each source pop- 
ulation in the above-defined region of the H.E.S.S. Galactic Plane 
Survey is shown in Table[5] 

6. Conclusion 

In this work, we search f or VHE counterpar ts of each Galactic GeV 
source in the OFGL catalog Abdo et al. ( 2009a), based on spatial coinci- 
dence. This study benefits significantly from the increased LAT angular 
resolution and its better sensitivity over previous instruments. 

Compared to the EGRET era, not only are there more coincident 
sources (improvement in quantity), but improvements in quality start to 
emerge. With the much better sensitivity of LAT, weaker sources are de- 
tected which were unknown in the EGRET era. New relations between 
the GeV and TeV spectra are revealed. A single spectral component is 
unable to describe some sources detected at both GeV and TeV ener- 
gies. Two spectral components may be needed in these cases to accom- 
modate the SEDs, where the VHE flux is higher than a PL extrapolation 
from GeV energies. 

A high fraction of Fermi bright sources are found to be spatially co- 
incident with a VHE y-ray source. This shows that there exists a com- 
mon GeV/TeV s ource population, a conclusion that is in stark disagree- 
ment with that of lFunk et aU d2008l) . 

Acknowledgements. We thank Emma Ona-Wilhelmi and Werner Hofmann for 
useful discussion. 
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Fig. 1. MeV to TeV spectra of four EGRET pulsars and their proposed nebulae. Upper left: Crab 
(OFGL J05 34.6+2201). The pul sar (dotted line) and nebula (dashe d-dotted line) spectral c omponents are those 
reported in lGrondin et al.1 (l2009h . T he VHE spectra ar e taken from lAharonian et al.l d200 6d). and the MILAGRO 
measurement at 20 TeV is shown (lAbdo et all 120071). Upper ri ght: Geminga (OFGL J0634.0+1745). The pul- 
sa r spectrum (dotted lin e) is that reported in ICelik et all ( I2009I) . The triangle denotes the upper limit r eported 
in lFinnegan et all (120091) . and the MILAGRO measurement at 20 TeV is also indicated (Abd o et alll2007l) Lowe r 
left: Vela (OFGL J0835.4-4510). The dotted line represents the Vela spectrum as shown in lAbdo et all (l2009d) . 
while the nebula comp onent is constrained by t he two triangles joined by the dashed-dotted line. The curved VHE 
spectrum is taken from lAharonian et al.l(l2006bl) . Lower ri ght: PSR B1706-44 ( OFGL J1709.7-4428). The dashed- 
dotted line denotes the two PL model spectrum derived in Garg ano et all (12009). Both LAT energy spec tra (though 
differe nt above 3 GeV) are consistent w ith the photon flux in the 1-100 GeV band of this source (Abdo et al., 
2009a). The VHE spectrum is taken fromHopp e et all (120091) . 
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Fig. 2. MeV to TeV spectra of the four radio pulsars first detected in y-rays by LAT. The spectra below 3 00 GeV 
are taken from Abdo et al] d2010t) . Up per left: The VHE spectr um is taken from lAharonian et ail d2006gl) . Upper 
right: The VHE spectrum presented in lAharonian et al.l d200 7b) is shown. The two curves represent the upper and 
lower limits of the s pectrum, taking measurement errors into account . Lower left : The fl ux at 35 TeV is taken from 
I Abdo et al.1 d2009gt). Lower righ t: The flux at 20 TeV is taken from lAbdo et alj d2007l) and the VHE spectrum is 
taken from lAcciari et all d2009d) . 
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Fig. 3. MeV to TeV spectra of the three new pu lsars discov e red in a blind search for pulsation in the LAT 
data. The spectra below 300 GeV are taken from lAbdo et al.l (l2010l) . Up per left: The VHE spectrum is taken 
fromlAharonian ~etal1 (l200 6g). Upper right: The VHE spectru m presented inlDiannati-Atai et al.l (l2008al) is shown, 
together with the c oincident MILA GRO source flux at 20 TeV (lAbdo et al.ll2007l) . Lower left: The VHE spectrum 
is that presented in lAharonian et al.l d2005bl) . while the MILAGRO flux at 20 TeV is also shown (lAbdo et aUl2007l) . 
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Fig. 4. MeV to TeV spec tra of four SNR/PWN ca ndidate OFGL sources. Spectra at >100 GeV e nergies are taken 
from lAlbert et all (l2007h (MAG IC J0616+225) lAharonian et all (120061 (HESS J1616-5081. [Aharonian et a"D 
(2008$ (HESS J1714-385), and lAharonian et al.1 (l2008d) (HESS J18 01-2331. The broken PL spectrum (dashed- 
dotted line) derived for 3EG J0617+ 2238 is taken fromlBertsch et alJ (120001) . The flux at 35 TeV at the position of 
OFGL J0617.4+2234 is that given in lAbdo et al.1 d2009gh . 
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| 0FGLJ1 834.4-0841 | 
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Fig. 5. MeV to TeV sp ectra of three SNR/PWN candidate 0FGL sources. Spectra at >100 GeV energies are taken 
from lAharonian et alj (120061 (HESS J1813-178 and HESS J 1834-0 87). For HESS J1923 +141, an assumed pho- 
ton index of 2. 8 is used in deriving the spectrum using the flux given in lFiasson et al. | (120091). and the flux at 35 TeV 
is that given in lAbdo et alJ (l200 9g). There is evidence of a steepening above several GeV (lAbdo et al.L l2009e). 
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Fig. 6. MeV to TeV spe ctra of four unident ified OFGL sources. Sp ectra at >100 GeV energies are taken from 
Ahar onian et all J2007ah (HESS J1023-5751. lAharonian et all (120061 (HESS J1634-472 and HESS J1804-216), 
and lTibollaetal](l2009bh (HESS J1741-302). 
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Fig. 7. MeV to TeV spe ctra of three unidentif ied OFGL source s . Spectr a at >100 GeV energies ar e taken from 
Aharo nianetail (120065) (HESS J 1837-069). lAharonian etail (l2008dl) (HESS J1841-055), and IChaves etaT 
(2009b) (HESS J1848-018). The flux at 35 TeV at the position of OFGL J1848.6-0138 is that given in lAbdo et al 
(20093. 
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Fig. 8. MeV to TeV spectra of the two X-ray binaries in the OFGL catalog. The phase-a veraged exponential cut-off 
spectr um in the GeV ra nge of LS I +61 303 {left) is taken from I Abdo eTal] J2009dl). Tha t of LS 5039 is taken 
from lAbdo et"ai] (l2009fl) . Spectra at >100 GeV energies ar e taken fr om Albert et alj (120091) (for a partial phase of 
LS I +61 303 during which VHE emission is detect ed) and Aharonian et alj d2006eh (for two phases of LS 5039). 
The flux at 35 TeV for LS I +61 303 is that given in lAbdo et alJd2009gl) . 
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OFGLJ1746.0-2900 | 




E(eV) 

Fig. 9. MeV to TeV spectrum of the Galactic Center. The VHE spectrum is taken from lAharonian e t al.l (T2 006h) 
while the broken PL spectrum (dashed-dotted line) derived for 3EG J1746-2851 is taken from iBertsch et all 
(2000). 




Fig. 10. Photon flux in the 100 MeV - 1 GeV band versus that in the 1-10 TeV band for 
coincident GeV/TeV sources. 
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Fig. 11. Photon flux in the 1-100 GeV band versus that in the 1-10 TeV band for coincident 
GeV/TeV sources. 




Photon index in the 1-10 TeV band 



Fig. 12. Photon index in the 0.1-100 GeV band (derived according to Sect. 13.31 ) 
versus photon index in the 1-10 TeV band, for coincident GeV/TeV sources. 
0FGL J1923.0+1411/HESS J1923+141 is not included since its VHE photon index is not 
known. 
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Notes. The top panel shows likely members of Pismis 11. The second panel contains likely members of Alicante 5. The bottom panel displays 
stars outside the clusters. 

<D) based on timing information {b) Source class according to Abdo et al. (2009a). PSR: pulsar; SNR/PWN: supernova remnant/PWN; HMXB: 
high-mass X-ray binary; Unid: unidentified sources. The classification of those sources as SNR/PWN is based on spatial coincidence only. 
(r) 95% positional error (rf) based on spatial coincidence only (e ' An entry of "PS" indicates that the source is point-like with respect to the point 
spread function of the respective instrument. u) These recent source discoveries are preliminary and they have been published in the referenced 
conference proceedings only. (s) The VHE emission has been detected by H.E.S.S. towards the direction of the massive stellar cluster. The 
coordinates refer to the nominal position of Westerlund 1. (/,) The extent of the source is not clear. The given value is est i mated from the radial 
profile shown in Fig. 4 of the corresponding reference. (,) The association is based on a morphological study (Ab do et al. , 20093). U) 

The extent 

of the source is not clear. The given value is estimated from the sky excess map in the corresponding reference. 



Re ferences. (1) lAharonian et alj d2006dh: ( 2) lAharonian et al.1 J2006bh: (3) lAharonian efaTJ j2006d) ; ( 4) iHoppe et all (2009) 
(5) lAharonian et alJ d2007bli: (6)lAharonian et alJ<2009h : ( 7)lAharonian et alj d2005bl) : (8 )lAcciari et all d2009bh : (9) lAharonian et al.l l2006f) 
(10) 
(15) 



Ohm et alj d2009h: (1 l)[Aharonian et alj d2008bl): (12) lAharonian et all teOOScK T (1 3)lFiasson et ^aU d2009l) ; (1 4) lAharonian et alJ (1 2007a) 
Aharonian et al. (2008d); (16) HooDe et al. (2008); (17) Chaves et al. (2009b); (18) Acciari et al. (2009a); (19) Aharonian et al. (2006e). 



Table 2. OFGL sources with coincident MILAGRO source, but without plausible coincident reported VHE y-ray sources. See Table[T]for the nomenclature. 
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Table 3. OFGL sources with a possibly coincident VHE y-ray source. See TableQ]for the nomenclature. 
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Notes. <a) These recent source discoveries are preliminary and they have been published in the referenced conference proceedings only. 

Refere nces. (1) lAharonian etal] d2006fl) ; (2) iTibolla et al.1 d2009bh ; (3) lAcero et al.1 j2O10h : (4) IWeinstein et"aH i |2009l) : (5) lAcciari et al] 
(2009cj). 



Table 4. H.E.S.S. sources which have a coincident OFGL source but do not have a 3EG counterpart as in Reim er et al. I (120081) 



H.E.S.S. sources OFGL sources coincidence level 
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Table 5. Number of coincidence cases for each source population (excluding extragalactic sources) in the region I = -85° to 60°, b = -3° to 3°. 



LAT Source class 


OFGL sources 


spatially 
coincident cases" 


pulsars 


10 


4 


SNR/PWN candidates 


11 


6(7) 


Unidentified sources 


19 


5(9) 


Total 6 


41 


16(21) 



Notes. (a) The numbers in brackets include possibly coincident cases (P). <6) including LS 5039 



